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ABSTRACT 


In the centers of dense star clusters, close encounters between stars and compact objects are likely to occur. We study tidal disruption 
events of main-sequence (MS) stars by stellar-mass black holes (termed uTDEs), which can shed light on the processes occurring 
in these clusters, including being an avenue in the mass growth of stellar-mass BHs. Using the moving-mesh hydrodynamics code 
AREPO, we perform a suite of 58 hydrodynamics simulations of partial / TDEs of realistic, MESA-generated MS stars by varying the 
initial mass of the star (0.5 Mo and 1 Mo), the age of the star (zero-age, middle-age and terminal-age), the mass of the black hole 
(10 Ma and 40 Mo) and the impact parameter (yielding almost no mass loss to full disruption). We then examine the dependence of 
the masses, spins, and orbital parameters of the partially disrupted remnant on the initial encounter parameters. We find that the mass 
lost from a star decreases roughly exponentially with increasing distance of approach and that a 1 Mo star loses lesser mass than a 
0.5 Me. Moreover, a more evolved star is less susceptible to mass loss. Tidal torques at the closest approach spin up the remnant very 
close to break-up velocity when the impact parameter is low. The remnant star can be bound (eccentric) or unbound (hyperbolic) to 
the black hole: hyperbolic orbits occur when the star's central density concentration is relatively low and the black hole-star mass ratio 
is high, which is the case for the disruption of a 0.5 Mo star. Finally, we provide best-fit analytical formulae for the aforementioned 


1. Introduction 
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> Tidal disruption events (TDEs) of stars by supermassive black 
CN holes (SMBHs) have been a subject of significant interest in the 
W) past decade (see Stone et al. 2019; Gezari 2021 for reviews). 
OO TDEs by SMBHs are observed as transients in multiple wave- 
lengths (Gezari 2021), with ~ 100 such events having been de- 
= tected by observatories in optical (ASAS-SN: Kochanek et al. 
CO 2017, ZTF: Bellm et al. 2019, PTF: Law et al. 2009, PS: Cham- 
m bers et al. 2016, ATLAS: Tonry et al. 2018, SDSS: Kollmeier 
Cf) et al. 2019, OGLE: Udalski et al. 2015), X-ray (ROSAT: Voges 
N et al. 1999, Swift: Gehrels et al. 2004, XMM: Jansen et al. 2001, 
> Chandra: Weisskopf et al. 2000) and UV (GALEX: Martin et al. 
° = 2005). 
In this work, we examine TDEs of stars by stellar-mass black 
d holes or SBHs (termed uU TDEs). These are less studied and have 
only garnered interest recently. Although uTDEs have lower ob- 
servable rates (Perets et al. 2016), they can shed light on the 
processes occurring in the centers of globular and nuclear star 
clusters where they are most likely to occur. In particular, they 
are an important avenue in the mass growth of SBHs to form 
intermediate-mass black holes (IMBHS, e.g., Stone et al. 2017; 
Rizzuto et al. 2023). However, these studies make simplistic as- 
sumptions about the mass accreted onto a black hole after a TDE, 
which highlights the need for detailed simulations to model this 
interaction more accurately for the next generation of globular 
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range of parameters that can be incorporated into cluster codes to model star-black hole interaction more accurately. 


Key words. stars: kinematics and dynamics — black hole physics — hydrodynamics — gravitation 


cluster simulations. Partial tidal disruptions (PTDEs), which are 
more likely than full tidal disruptions (FTDES), can be responsi- 
ble for the tidal capture, and subsequent encounters, of the rem- 
nant star by the BH if the remnant ends up in a bound orbit. The 
remnant stars of such interactions tend to be spun up from the 
torque due to the BH and can have peculiar internal structures 
(e.g., Alexander & Kumar 2001). 


HTDEs are also of interest in the context of low-mass X- 
ray binaries (LMXBs) and the newly-discovered Gaia BHs (El- 
Badry et al. 2023a,b; Chakrabarti et al. 2023). Isolated binary 
evolution cannot satisfactorily explain the observed rates of bi- 
naries with BHs and low-mass stellar companions (e.g., Podsi- 
adlowski et al. 2003; Kiel & Hurley 2006). Hence, the dynam- 
ical formation of such binaries in clusters, through PTDEs and 
subsequent tidal captures, may very well be a crucial process to 
explain them. 


Accurate modeling of post-disruption orbits and the internal 
structure of remnant stars requires detailed hydrodynamics sim- 
ulations. The first detailed study on y TDEs was carried out by 
Perets et al. (2016). They estimated that these events could occur 
in Milky Way globular clusters at a rate of 105 yr! MWGal !, 
and might observationally resemble ultra-long GRB events. Sim- 
ilar rates were obtained for scenarios where a supernova natal 
kick to a newly formed BH results in a chance encounter with its 
binary companion. Wang et al. (2021) used smoothed-particle 
hydrodynamics (SPH) to find the fallback rates onto the BH 
in the case of partial /TDEs of polytropic stars. Kremer et al. 
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(2022) performed a large suite of SPH simulations of partial- 
and full-disruptions, with varying BH and stellar masses, stellar 
polytropic indices, and impact parameters. Among other things, 
they observed that stellar structure plays a crucial role in the 
boundedness of the remnant stars after partial disruption, with 
the possibility of less centrally ‘concentrated’ stars ending up in 
hyperbolic orbits. Xin et al. (2024) looked at low-eccentricity 
‘tidal-peeling’ events of realistic MS stars, with the star being 
slowly stripped away in multiple orbits. 

Hydrodynamics studies on three-body encounters involving 
SBHs in globular clusters have also garnered interest, motivated 
by the significantly larger cross-sections of binary-single tidal 
interactions. Lopez et al. (2019) studied the interaction of binary 
BHs and polytropic stars and showed that tidal disruption can 
alter the spin of one of the BHs. An extensive series of studies 
was carried out by Ryu et al. (2022, 2023a,b, 2024b), using both 
SPH and moving-mesh codes, on the different combinations of 
close encounters between realistic MS stars and SBHs. 

In our study, we use the moving-mesh code AREPO (Springel 
2010; Pakmor et al. 2016; Weinberger et al. 2020) to simu- 
late uTDEs of low-mass MS stars with SBHs. It should be 
noted that Ryu et al. (2023a,b, 2024b) also employed AREPO 
to model TDEs with remarkable success. We generate realis- 
tic stellar models from detailed 1D MESA (Paxton et al. 2011) 
profiles. TDEs of MESA-generated stars by SMBHs and IMBHs 
have been investigated in the past (e.g., Gallegos-Garcia et al. 
2018; Golightly et al. 2019; Law-Smith et al. 2019, 2020; Ryu 
et al. 2020a,b,c,d; Kıroğlu et al. 2023), but uyTDE studies have 
typically probed polytropic stars. In this paper, we focus on the 
dependence of post-disruption mass, spin, and orbital parame- 
ters on the initial conditions and provide best-fit functions for 
the same. These fits can be incorporated into N-body or other 
cluster codes for better treatment of star-BH encounters. 

The paper is organized as follows. In Section 2, we detail our 
simulation methodology, including our grid of initial conditions. 
We briefly describe the analysis of certain quantities in Section 
3 and present our results and best-fit functions in Section 4. Sec- 
tion 5 and Section 6 discuss the implications of our work and 
conclude, respectively. 


2. Methods 
2.1. Hydrodynamics 


We simulate the disruption of stars by black holes using the 
moving-mesh magnetohydrodynamics code AREPO (Springel 
2010; Pakmor et al. 2016; Weinberger et al. 2020). AREPO is a 
massively parallel 3D magnetohydrodynamics code with gravity 
that inherits many advantages of the two popular schemes, Eu- 
lerian grid-based finite-volume codes and Lagrangian smoothed 
particle hydrodynamics (SPH). Although initially developed for 
cosmological simulations, AREPO has been successfully em- 
ployed in phenomena involving stars, e.g., tidal disruptions and 
encounters (Goicovic et al. 2019; Ryu et al. 2023a,b, 2024b), 
TDEs in active galactic nuclei disks (Ryu et al. 2024c), collisions 
of main-sequence stars (Schneider et al. 2019), giant stars (Ryu 
et al. 2024a), white dwarfs (Pakmor et al. 2013, 2021, 2022; 
Burmester et al. 2023; Gronow et al. 2021; Glanz et al. 2023) and 
neutron stars (Lioutas et al. 2024), and common envelope evo- 
lution (Ohlmann et al. 2016a,b; Sand et al. 2020; Kramer et al. 
2020; Ondratschek et al. 2022). 

To solve the fluid equations, AREPO builds an unstructured 
Voronoi mesh with cells of varying volume and calculates fluxes 
between the cells in a finite-volume approach. The Voronoi mesh 
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Table 1: Parameters of the Z, = 0.02 MS star models. 


my [Mo] ts [ Gyr] H, ry [Ro] ldyn,x [ hr] Prone 
0.5 0.01 0.70 0.44 0.19 0.46 
0.01 0.70 0.90 0.38 0.29 

1.0 4.68 0.34 1.02 0.46 0.20 
8.60 0.00 1.22 0.60 0.12 


Notes. The columns indicate stellar parameters - the mass m,, the age 
tx, the central hydrogen abundance H,, the radius r,, the dynamical time 
fayn,» and the inverse of the density concentration parameter p... 


moves in time according to the approximate bulk velocities of 
the fluid elements. Gravity is handled in a tree-particle-mesh 
scheme (see Bagla 2002), with the minimum gravitational soft- 
ening for the gas (star) cells set as one-tenth of the smallest gas 
cell for accuracy. AREPO allows for adaptation of spatial resolu- 
tion according to arbitrary criteria on top of the default adaptiv- 
ity to density, inherited from the near-Lagrangian nature of the 
scheme. It also includes particles that interact only gravitation- 
ally. We assume Newtonian gravity and do not include magnetic 
fields in our simulations. 


2.2. Stellar models 


We use the 1D stellar modeling code MESA (Paxton et al. 2011) to 
generate main-sequence (MS) stellar models with initial masses 
m, = 0.5 Mo and 1.0 Mo, and metallicity Z, = 0.02 (near-solar) 
at different ages. These masses are typical for globular clusters 
that predominantly consist of old, low-mass stars (e.g., Salaris 
& Weiss 2002; De Angeli et al. 2005; Marin-Franch et al. 2009), 
while the metallicity is higher than those of most clusters, which 
typically have a tenth (or less) of Solar metallicity (e.g., Harris 
2010)). However, the effect of metallicity is not expected to be 
significant and is quantified in Section 4. MESA solves the stel- 
lar structure equations to compute the evolution and provides us 
with stellar mass- and radial-profile parameters (e.g., densities, 
temperatures, and chemical abundances) at different evolution- 
ary phases. We convert these 1D profiles to 3D AREPO initial 
conditions using the scheme adopted by Ohlmann et al. (2017), 
with the Helmholtz equation of state (Timmes & Swesty 2000). 
After generating the star, we relax it for five stellar dynamical 
timescales ftaynx = (r /Gm,)0?, where r, is the stellar radius 
and G is the gravitational constant. We then use this relaxed star 
for the disruption simulations. 


Stellar density profiles are a major factor in determining 
whether a star undergoes a partial tidal disruption (PTDE) or a 
full tidal disruption (FTDE) for a given periapsis distance (e.g., 
Ryu et al. 2020c). To that end, we consider MESA models of 
three evolutionary stages of 1.0 Mo MS stars with varying core 
hydrogen abundances H,, which correspond to distinct density 
profiles (see also Gallegos-Garcia et al. 2018; Golightly et al. 
2019; Law-Smith et al. 2019; Goicovic et al. 2019). The first is 
close to the onset of hydrogen burning (He ~ 0.70), i.e., zero-age 
main-sequence (hereafter ZAMS), the second is approximately 
midway through the main-sequence (H, ~ 0.34), i.e., middle- 
age main-sequence (hereafter MAMS), and the third is close to 
the depletion of core hydrogen (He « 0.00), i.e., terminal-age 
main-sequence (hereafter TAMS). The total stellar mass remains 
nearly constant throughout the MS lifetime owing to the insignif- 
icant wind mass loss. In the case of the 0.5 Mo stars, we exam- 
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Fig. 1: Density (left) and temperature (right) profiles of a 1 Mo MS star of different models (see Table 1). The dashed and solid 
lines represent the 1D profiles generated by MESA and initialized (and subsequently relaxed for five dynamical timescales) in 3D in 
AREPO respectively. They agree very well over most of the range except close to the stellar surface because of smoothing during 
relaxation in AREPO. An older MS star (with lower H,) has a denser and hotter core and a puffier outer layer than a younger MS star 
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Fig. 2: Similar profiles as Figure 1, but for a 0.5 Mọ MS star. A 0.5 Mo star of age ~ 13.5 Gyr with H. = 0.60 (purple curves) is 
very similar to a ZAMS star with H, = 0.70 (red curves), albeit with a slightly larger radius. 


ine only ZAMS profiles (H. z« 0.70). This choice is motivated by 
the fact that their internal structure barely changes over time (see 
Figure 2) owing to their large MS lifetimes. One could choose 
a 0.5 Mo star at an age comparable to a Hubble time (e.g., Ryu 
et al. 2020a,c) but the results are expected to be very similar. Ta- 
ble 1 lists the stellar parameters of the MS star models for the 
stellar masses and ages we simulate. 


Furthermore, we performed resolution scaling tests to ensure 
that the PTDE simulations were robust. To this end, we varied 
the number of fluid (star) cells making up a 1 Mo star — 5 x 104, 
8 x 10+, 2 x 105, 5 x 105, and 8 x 10°. We then performed PTDE 
simulations involving a 1 Mo star and a 10 Mo BH for each res- 
olution. The masses lost due to tidal disruption plateaued close 
to the resolutions 2 x 10? and 5 x 10°, with the mass loss val- 


ues using 5 x 10? and 8 x 10? cells being essentially the same. 
Hence, we chose an initial resolution of 5 x 10? fluid cells for 
the stars in this study. It should be noted that the number of cells 
increases as the simulations progress since AREPO automatically 
adds extra cells when the density gradient is high. 


Figure | shows that more evolved 1 Mọ MS stars have a 
denser (left panel) and hotter (right panel) core and a puffier 
outer layer than less evolved MS stars. Since the three stars have 
different radii rą, they also have different fayn,» and different tidal 
radii r, = (mgy/m,)! Pr, for the same BH mass mpy. Figure 2 
shows that 0.5 Mo stars of different ages - ZAMS (H. = 0.70) 
and ~ 13.5 Gyr (He = 0.60) — have very similar density and tem- 
perature profiles, which justifies excluding the latter from our 
simulation runs. It should be noted that the AREPO profiles, espe- 
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Table 2: Initial parameters of our suite of simulations. 


mx [Mo] mgu [Mo] He b= rp/ 
0.5 10.0, 40.0 0.70 0.25, 0.33, 0.50, 0.75, 1.00, 1.50, 2.00, 2.50 
1.0 10.0, 40.0 — 0.70, 0.34, 0.00 0.25, 0.33, 0.50, 0.75, 1.00, 1.50, 2.00 


Notes. The columns indicate the initial parameters - the stellar mass m, , the BH mass mpy, the core hydrogen fraction H., and the impact parameter 
b. In each row, every comma-separated combination of values in all columns is simulated, to give a total of 58 simulations. 


cially the temperature, diverge significantly from the MESA ones 
near the surface of the star (corresponding to a fractional mass of 
~ 1077-1072). This amount of mass near the surface that devi- 
ates from the MESA model sets the minimum amount of stripped 
mass in PTDEs that we can resolve. Therefore, we only present 
the results of PTDEs yielding a mass loss = 10^^m, in this pa- 
per. 


2.3. Black hole 


The black holes are initialized in AREPO as non-rotating sink par- 
ticles that interact gravitationally and grow in mass through ac- 
cretion. We set the gravitational softening length of the BH to be 
ten times the minimum softening length of the gas (star) cells. 
The scheme we use for accretion onto the BH, described here 
briefly, is the same as the one used in Ryu et al. (2022, 2023a,b, 
2024b). Firstly, the cells around the BH within a radius of 10frg, 
where gravitational radius rą = Gmpy/c’, are identified. Sec- 
ondly, the accretion onto the BH is calculated using a weighted 
average of radial flux to account for the higher contribution of 
gas closer to the BH and vice-versa. Thirdly, mass is extracted 
from the neighboring cells of the BH, and momentum is inserted 
in the BH to conserve these quantities. Finally, the BH is spun up 
from accretion by adding angular momentum from the selected 
neighboring cells. The radiation feedback from accretion is not 
taken into account. Although angular momentum is conserved 
through this scheme, there is no guarantee that the total energy 
is conserved. However, given that accretion onto the black hole 
is small in our simulations, the energy lost due to accretion is 
negligible. 

We run simulations with SBH masses of 10 Mo and 40 Mo, 
which are in proximity to the lower- and higher-mass peaks of 
observed LIGO BH masses (e.g., Abbott et al. 2021, 2023). 
Given the mass of the BH mgg, the mass of the star m,, and 
the radius of the star r,, the classical tidal radius is defined as 
r, = (nga /m,) Pra. 


2.4. Initial conditions 


The distance of the closest approach is characterized by the im- 
pact parameter: b = ry/r ' where ry is the periapsis distance of 
the initially parabolic orbit. We vary the mass of the star: 0.5 Mo 
and 1 Ms, and the mass of the black hole: 10 Mo and 40 Mo. 
We then generate a suite of simulations of b varying from 0.25 
to 2.00 (1 Mo stars) or 2.50 (0.5 Mo stars), and the three stellar 
ages (density profiles) to determine the post-disruption orbital, 
mass and spin parameter. The additional simulations for 0.5 Mo 
and b = 2.50 are performed to properly analyze the PTDE trends 
since 0.5 Mo stars are fully disrupted for b x 1.00. 


! B = b! is often introduced in the literature for TDEs by SMBHs as 
the ‘penetration factor’. 
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We also assume that the initial stars are non-rotating and that 
the star-BH encounters are parabolic. This can be justified for 
encounters in globular clusters because the velocity dispersions 
of Milky Way globular clusters are in the range ~ 1-20 km s~! 
(e.g., Baumgardt & Hilker 2018). For instance, given a velocity 
dispersion o = 10km sl, |I- e| = 10-4 (very close to parabolic) 
for the parameters considered in this paper. 

Each simulation is run for ~ 70 tayn,» after periapsis passage. 
This duration was chosen to ensure that the stellar masses reach 
steady values for proper analysis. Table 2 provides an overview 
of the parameters we use in our suite of simulations. 


3. Analysis 
3.1. Calculation of bound and unbound mass 


To find the center of mass of a self-bound object (i.e., star before 
disruption, remnant after disruption) in each snapshot, we used 
an iterative procedure fairly similar to Guillochon & Ramirez- 
Ruiz (2013), with some differences to improve the accuracy of 
the identification of bound and unbound cells. For each snapshot, 
we chose the gas cell with the maximum density as the ‘initial’ 
guess for the center of mass (CoM) position. Subsequently, we 
calculate the specific total energies of all the star cells relative 
to the star CoM, &cen,com, and the BH, &cen.pH, as the sum of 
their relative kinetic, potential and internal energies. We consider 
a cell to be bound to a self-gravitating object if Ecen,com < 0, 
Ecell,CoM < Ecell,BH- The last condition ensures that the cells bound 
to the star are ‘close’ to the CoM. We compute a ‘new’ CoM 
position (and velocity) using these bound cells and determine 
the ‘new’ specific energy. We iterate this process until there is 
no change in the relative position of the CoM. Subsequently, we 
calculate its ‘final’ bound mass. In the case of an FTDE, we also 
visually inspect the surface density plots and assign the mass 
bound to the star to be zero. 

The mass bound to the BH is calculated in a similar fash- 
ion, though the position of the BH is trivially known from the 
simulation. Finally, we consider any cells with €cen,com 2 0 and 
Ecell,Com = 0 unbound from the system. 


3.2. Calculation of orbital and spin parameters 


For each snapshot post-disruption, we compute the instanta- 
neous Keplerian orbital parameters — the semimajor axis a and 
orbital eccentricity e — using the current values of the mass 
bound to the star m,, the mass of the BH (plus the gas mass 
bound to it) mgu, the positions and velocities of the star’s CoM 
(com; Vcom) and the BH (rgy, Vgy). For completeness, with 
r = com — fg and v = Vcom — Vau, the equations are as follows: 


2 
a=( y (1) 


-1 
r G(m, + 21 
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Table 3: Evaluated values of the impact parameter for the full 
disruption of a star with uniform density, beroe = Derpg/Pconc- 


y  brrpE Pohe PeIDE 
4/3 0.50 0.26 1.92 
5/3 1.08 0.55 1.96 


Notes. The columns indicate the polytropic index y, the impact param- 
eter for full disruption bfrpe, the inverse of the density concentration 
factor PES and the aforementioned parameter bgrpg. The bgrpg values 


for polytropic stars are from Mainetti et al. (2017). 


zit 


We determine the spin angular momentum L, about the 
star's CoM using the masses m;, relative positions r; and veloci- 
ties v; of the bound gas cells. In the cases of full disruption, we 
ignore the orbital and spin parameters. 

The ‘final’ post-disruption mass, orbital, and spin parameters 
reported in the following section are the means and standard de- 
viations of these quantities during the last ten snapshots of each 
simulation (the time between consecutive snapshots was chosen 
to be approximately equal to the initial dynamical timescale of 
the star). In the cases where a partially disrupted star returns on a 
second passage within the simulation time, we chose snapshots 
close to the apoapsis of the first passage for the parameter calcu- 
lations. 


(r x y? 0.5 
G(m, + mpy)a 


(2) 


4. Results 


Figure 3 shows selected snapshots of our suite of simulations. It 
can be visually seen that the final outcomes (FTDE or PTDE) 
depend on all the parameters that we considered — the stellar 
mass m,, the core hydrogen fraction He, the BH mass rpg, and 
the impact parameter b. 

The following sections detail the quantitative results. It 
should be noted that, henceforth, initial star parameters are given 
the subscript x (e.g. m,, L,), while the post-disruption remnant 
parameters are given the subscript o (e.g. ms, Lo). Table A.1 pro- 
vides detailed initial and post-disruption parameters for our 58 
simulations. 


4.1. Effect of density concentration factor 


Ryu et al. (2020b) showed that, in the case of SMBHS, the criti- 
cal impact parameter for full disruption, bprpe, is inversely pro- 
portional to the density concentration factor Peone = (p. /p)!^, 
where p and p, are average and central stellar density respec- 
tively (see also Law-Smith et al. 2020). A larger value of (conc 
represents a denser core with a puffier outer layer, and vice- 
versa. This indicates that a star with higher core density requires 
a smaller impact parameter to fully disrupt. Quantitatively, berpg 
can be written as: 


(3) 


Here, bprpe is the critical impact parameter for the full dis- 
ruption of a hypothetical star with uniform density, i.e., Oconc = 1 
(see also Equation 15 of Ryu et al. 2020b). Table 3 displays the 
values of bfrpe of stars with polytropic indices y (where P œ p”) 


cum T 
DeTDE = Poe PFTDE 


4/3 and 5/3 as calculated by Mainetti et al. (2017). Shown also 
are the analytically computed (conc for these polytropes, and the 
subsequently evaluated (using Equation 3) bprpe values. We also 
estimate bprpg to be 1.95 + 0.20 from our suite of simulations 
by fitting the remnant masses for given initial stellar and BH 
masses (see the following section for details), and find the val- 
ues of bprpg. Despite the significant difference in the masses 
of the black holes, our estimated values agree with those from 
Mainetti et al. (2017). Estimates of this factor calculated by Ryu 
et al. (2020b), Law-Smith et al. 2020 and Kıroğlu et al. (2023) 
all lie within the error range. Since the former two studies are 
on SMBHs and the latter is on IMBHs, we reiterate that bfprpe 
is almost independent of the mass of the BH for a very large 
range of BH masses. The consensus between the previous works 
in spite of different numerical schemes (general relativistic grid- 
based hydrodynamics, adaptive-mesh refinement, and SPH re- 
spectively) and the current work strengthens this point. There- 
fore, we will use our computed value of bprpg = 1.95 + 0.20 for 
the fit functions in the following sections. 


For completeness, Figure 4 (and Table 1) shows the trend 
pa. for MESA 1 Mọ MS stars, with our three stellar models 
highlighted. We see that a more evolved star has a smaller p... 
(harder to fully disrupt) compared to a less evolved star (eas- 
ier to fully disrupt) and that this relation is close to linear as a 
function of stellar age tą. Figure 5 shows the averaged o1. val- 
ues for a range of masses of MESA-generated MS stars (0.1 Mo 
— 20.0 Mo) of two metallicities — 0.02 (near-Solar) and 0.002 
(sub-Solar). We see that higher mass MS stars have denser cores 
and puffier outer layers (hence lower p; values, which plateau 
for stars of masses z 1 Mo at ~ 1.6-2.0), while the opposite is 
true for lower mass MS stars. Moreover, metallicity affects pz), 
only slightly, with the metal-poor stars being slightly less (more) 
centrally concentrated when m, > 1.0 Mo (m, x 1.0 Mo). The 
(approximate) values of pl. are crucial to applying our fits (see 
the following sections). 


4.2. Post-disruption stellar and BH masses 


After a disruption event, a star loses mass, some of which gets 
bound to (and eventually accreted onto) the BH while the rest of 
it becomes unbound from the system. Evidently, the closer the 
approach of the star is to the BH, the larger the mass lost from 
the star. 


Figure 6 shows this trend for the case of a 1 Mo (left) and a 
0.5 Mo (right) star that is disrupted by a 10 Mo (top) and a 40 Mo 
(bottom) BH. The post-disruption fractional mass loss Am, /m.,, 
where Am, is the difference between the star's initial mass m, 
and the remnant mass ms. It decreases almost exponentially with 
the impact parameter b. Moreover, a more evolved MS star loses 
less mass compared to a less evolved star for the same value 
of b. This is because a more evolved star, with a higher central 
concentration, requires a stronger tidal force (or smaller b) to 
strip the same amount of mass when compared to a less evolved 
star. Consequently, brrpg, corresponding to Am,/m, = 1, is also 
lower for more evolved stars. 


Furthermore, we see that roughly half the mass lost from the 
star stays bound to the BH, although this is not always true, espe- 
cially when there is little mass loss (see Table A.1). This bound 
mass decreases through time due to continuous debris interac- 
tion resulting in mass unbinding, and hence is an upper limit for 
the mass that can be accreted onto the BH. Details of the accre- 
tion process and feedback from the BH might prevent some of 
the mass bound to the BH from being accreted, but following 
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Fig. 3: Grid of log density slices of stars ~ 20 tayn after undergoing uTDEs, with the BH at the center of each panel. Each row is 
for a different mgu (top and bottom half for 10 Mo and 40 Mo respectively), m, or Hc, whereas the columns represent increasing b 
(from left to right). The dashed and solid circles around the BH, most visible in the rightmost panels, denote the tidal radius, rų, and 
the periapsis distance, rp, respectively. Note that the spatial range of each panel is not the same. 


this process is beyond the scope of our paper. The other half of Here, bprpg is as given in Equation 3. Since pcs. is taken 

the stripped mass is unbound from the BH-star system. into account, this fit differs from those obtained by Guillochon 

We can fit (with standard 1c errors) the fractional mass loss & Ramirez-Ruiz (2013) and Ryu et al. (2020c) for TDEs due 

Am, well by: to SMBHs. The remnant mass fit function does not explicitly 

Àm depend on m, and mg, but only on b and pco. Since possc 1s 

logio (=) = min {(—2.10 + 0.10) (b — berpg) , 0} (4) higher for more evolved stars, they have a smaller bprpg and 
* 
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Fig. 4: MESA-derived inverses of density concentration factors 
Daly, for a 1 Mo MS star as a function of stellar age. The three 
dashed vertical lines represent the ages of three chosen MS mod- 
els in our study, corresponding to the profiles in Figure 1 and the 
values in Table 1. The dependence with age is very close to lin- 
ear. 
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Fig. 5: MESA-derived inverses of density mean concentration fac- 
tors Pde for MS stars of a range of masses (0.1, 0.2, 0.5, 1.0, 2.0, 
5.0, 10.0 and 20.0 in units of Mo) and two metallicities — 0.02 
and 0.002. The shaded regions illustrate the variation in pil. 
over the MS lifetime, with a TAMS star having the lowest o... 
value and a ZAMS star having the highest ( the values are for 
ZAMS stars only in the cases of stars of masses « 1 Mo). The 
two dashed horizontal lines represent the corresponding values 
for polytropic stars, shown in Table 3. 


hence, less fractional mass loss. Figure 6 plots these best-fit 
curves, with the error bars shaded. It should be noted that this 
trend may not hold for larger values of b. 


4.3. Post-disruption stellar spins 


An initially non-rotating star gains spin after a tidal encounter 
due to tidal torques from the BH. Although we do not delve 
into the details of the spin-up, this results in differential, and 


not rigid, rotation of the star. It should be noted that, due to 
spurious motions in our initial AREPO stellar models, the spin 
angular momenta L, are not zero but ~ 1-2 x 107 gcm? s^! 
(~ 2x 1076 g cm? s7!) for 1 Mo (0.5 Mo) stars. Since these val- 
ues correspond to negligible equatorial surface velocities? of 
~ 10ms"!, the stars are considered to be non-rotating. The spin 
angular momentum of the remnant L, depends on the internal 
structure of the original star and the impact parameter. 

Figure 7 shows the trend of L, with respect to the impact 
parameter b for the case of a 1 Mo (left) and a 0.5 Mo (right) 
star being partially disrupted by a 10 Mg (top) and a 40 Mo 
(bottom) BH. In addition, we indicate the approximate, order- 
of-magnitude estimates of the break-up angular momenta (as- 
suming uniform density and rigid rotation) of the remnant stars, 
Lo break = 0.4 Moro Vo break 5. Here, we crudely estimate the rem- 
nant radius as (r,/ Ro) = (m,/ Ma" and assume the break- 
up velocity to be the Keplerian velocity at rx, Le., Vo break = 
vGm. [ r,. In reality, the remnant is significantly puffed up, and 
the break-up velocity would depend on the true moment of iner- 
tia of the remnant (and thus, its density profile) and the differen- 
tial rotation due to spin-up. For our range of simulation parame- 
ters, Lo remains below Ls break, although significant mass loss can 
bring the remnant very close to break-up, e.g., m, = 0.5 Mo and 
b € 1.00. For higher values of b (but still less than b ~ 2), L, de- 
creases almost exponentially as b increases. On the other hand, 
when b is less such that a fractional mass loss is > 2 x 107!, 
the trend reverses until FTDE. This trend reflects the counter- 
balance between mass loss and spin-up: for stronger PTDEs, al- 
though the remnant is spun up more rapidly, the remnant mass is 
smaller. 

The large spin-up for smaller b indicates that the specific an- 
gular momentum monotonically increases as b decreases. Mo- 
tivated by this, we considered the scaled (moment of inertia- 
adjusted) parameter Lo scaled = Lo/(Mo.frac T5 rao): Here, the frac- 
tional remnant mass rn; frac = Mo /M and the post-disruption es- 
timated fractional radius 7, frac/ Ro = mot e is a crude estimate 
of the radius of an MS star. 

We found the following exponential fit (with standard 1c er- 
rors) describes Lo scale Well: 


Le scale 
logo [A a ) = (-1.50 + 0.15) b + (52.00 + 0.30) (5) 
gcm? s 
From the definition of 7, frac before, we have: 
24 24 
Ms Am. 
L, = Lo scaled (=) = Ls scaled ( mi (6) 
my m, 


Here, (Am, /m,) is estimated using Equation 4. This equation 
implies that the post-disruption spin depends significantly on the 
fractional mass loss. Moreover, there is no dependence on mpy. 
Figure 7 plots these best-fit curves, with the error bars shaded. 


4.4. Post-disruption orbital parameters 


Depending on the stellar structure and the BH mass, a star in 
an initially parabolic orbit ends up in an eccentric or even a hy- 
perbolic orbit after a PTDE (see also Ryu et al. 2020c; Kremer 
et al. 2022; Kiroglu et al. 2023). The specific orbital energy of 


? Tn contrast, the Sun, a very slow rotator, has an equatorial surface 
velocity of ~ 2km s^!. 

? This is because hypothetical star (sphere) of uniform density has a 
moment of inertia of /, uniform = 0.4 m,r2 about an axis of rotation 
passing through the center. 
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Fig. 6: Post-disruption fractional mass loss, Am,/m,, for a star of mass 1 Mo (left) and 0.5 Mo (right), due to a BH of mass 10 Mo 
(top) and 40 Mo (bottom), as a function of impact parameter b and MS age (H. abundance). The dashed line at the top signifies the 
full disruption of the star. The fractional mass losses are largely independent of BH mass. Mass loss is higher for lower b values, and 
this decreases roughly exponentially (best-fit lines and standard error regions shown) with increasing b. A less evolved star (higher 
He) loses more mass for the same value of b. A 0.5 Mo star has a higher mass loss than a 1 Mo star for the same b. 


the star-BH system, £orb, is the sum of the specific energy in- 
jected into stellar tides, jac, the specific binding energy of the 
ejected material, €ping, and the specific ‘kick’ energy, Eick (see 
Kremer et al. 2022). Manukian et al. (2013) estimated the poten- 
tially unbinding ‘kick’ received by a star undergoing PTDE due 
to an SMBH from asymmetric mass loss and observed that it de- 
pends solely on the mass loss and not on the star-SMBH mass 
ratio. However, this is not true for uTDEs — the mass ratio de- 
termines the kick, and hence the boundedness, of the star-SBH 
system (see also Figure 2 of Kıroğlu et al. 2023 for TDEs due to 
IMBHs). 


Figure 8 shows the post-disruption specific orbital energy 
of the star-BH system, £o», for the disruption of a 1 Mo and a 
0.5 Mo star by a 10 Me (top panel) and a 40 Mo (bottom panel) 
BH. The values are normalized to the square of & = IlOkms"!, 
the typical velocity dispersion of Milky Way globular clusters 
(e.g., Baumgardt & Hilker 2018). More negative (positive) val- 
ues represent more bound (unbound) orbits. In general, £o be- 


Article number, page 8 of 17 


comes more negative with decreasing b. If the density concentra- 
tion factor Pconc is high enough (as in the case of a TAMS 1 Mo 
star; see Figure 4), this trend continues till FTDE. However, in 
the case of a ZAMS 1 Mo star (and less extremely in a MAMS 
star), there is a certain value of b at which £or» starts to increase 
with decreasing b. This is due to the momentum kick from asym- 
metric mass loss, which becomes significant when the mass loss 
is high. In the case of a PTDE of a 1 Mo ZAMS star by a 40 Mo 
BH, the post-disruption orbit can be hyperbolic (unbound) for b 
values close to FTDE. This trend is more extreme for a 0.5 Mo 
star, with hyperbolic orbits being possible even in the case of 
PTDE by a 10 Mo BH. Moreover, their velocities at infinity are 
high enough (~ 100—300 km s^!) to escape the cluster entirely. 
The two different BH masses result in quite different quantitative 
orbital parameters. Similar to Kremer et al. 2022, we can write 
each of the energy components as follows: 
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Fig. 7: Post-disruption spin angular momentum, L,, for an initially non-rotating star of mass 1 Mo (left) and 0.5 Mo (right), due to 
a BH of mass 10 Mo (top) and 40 Mo (bottom), as a function of impact parameter b and MS age (H, abundance). The dashed lines 
indicate the (approximate) estimated break-up angular momenta of the remnants Ls break respectively. The spin angular momenta 
are largely independent of BH mass. It decreases almost exponentially for sufficiently high b values, and it drops for very low 5 
values when the remnant mass is small (best-fit lines and standard error regions shown). The drop-off of angular momentum occurs 


at larger b for a less evolved star. 


Gm2 6 8 
Eide = —-—= (2) Ta + (=) T3 (7) 
Fyme Fp Tp 
GAm, 
£bind = 77 (8) 
Vy 
Gmggy Am; r? 
Ekik = —— 3—— (9) 
riis 


Here, T and T5 are terms which depend on the stellar mass, 
radius, and density structure, BH mass, and distance of approach 
(Press & Teukolsky 1977; Lee & Ostriker 1986). For the range 
of parameters considered, with mpg > m, and ry ~ rx, the val- 
ues of T» and T; lie in the range 0.01 — 0.1. These values are 
systematically lower for stars with lower pconc values. As a re- 
sult, in 0.5 Mo stars (and ZAMS 1 Mo stars), &ijge is dominated 
by €xick, resulting in them becoming unbound from the BH. It 


should be noted that eig. depends very strongly on the distance 
of approach. 


In addition, we find that the change in the specific orbital 
angular momentum of the star-BH system, hop, is not signif- 
icant. For orbits very close to being parabolic, which is the 
case for all our post-disruption orbits (see Figure 9), how « 
2am, +m) * [2Gmsur,. This implies that the post- 
disruption periapsis distance rp = a(1 — e), where a is the semi- 
major axis and e is the eccentricity, is nearly the same as the 
initial periapsis distance rp = br, which is also seen in PT- 
DEs by SMBHs (Ryu et al. 2020c). Consequently, given that 
a = —0.5G(m, + mpy)/Eorw, the calculation of e is straight- 
forward. Figure 9 shows the post-disruption e, computed using 
Equations | and 2, as a function of b which follows the trend in 
Eorb- 
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Fig. 8: Post-disruption normalized specific orbital energies of the star-BH system, £orb, for a star of mass 1 Mo (left) and 0.5 Mo 
(right) initially in a parabolic orbit, due to a BH of mass 10 Me (top) and 40 Mo (bottom), as a function of impact parameter b and 
MS age (H. abundance). The dashed line represents parabolic orbits and separates the parameter space into bound (eccentric) and 
unbound (hyperbolic) orbits. When mass loss is relatively low, a lower value of b generally results in more negative €o, in the case 
of 1 Mo star. This trend continues for the TAMS star. However, for the ZAMS and MAMS stars, significant mass losses can reverse 
this trend, with hyperbolic orbits also being a possibility (e.g., 1 Mo ZAMS star and 40 Mo BH mass with b < 1.00). 0.5 Mo stars, 


on the other hand, can become unbound for larger b, with boundedness also depending on the mass of the BH. 


4.5. Rates of tidal encounters 


The rate of uL TDEs per year per Milky Way-like galaxy (similar 
to Perets et al. 2016) can be expressed as: 
RutDE = NetusVx,c NBH, Lo y; (10) 
Here, Naus ~ 150MWGal ! is the number of globular clus- 
ters (in the Milky Way, e.g., Harris 2010), Nx, is the number 
of stars in a cluster core, Ngy, is the number of BHs in a clus- 
ter core, X = nrl + v? /0?) is the gravitational focused en- 
counter cross-section with periapsis distance ry = br; (see Sec- 
tion 2) and periapsis velocity vp = ~G(mpy + m,)/rp, and c 
is the velocity dispersion of the cluster. For globular clusters in 
the Milky Way Galaxy, the typical values of these parameters 
are Nyc = NycVz! ~ 107-107 pc? (e.g., Baumgardt & Hilker 
2018), Ngu, ~ 10-100 (e.g., Morscher et al. 2015; Askar et al. 
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2018), 7 ~ 1-20 km s^! (e.g., Baumgardt & Hilker 2018). With 
these values, and assuming vp >> c in globular clusters, and 
typical star and BH masses of 1 Mo and 10 Mo respectively, the 
rate is: 


= = = Notus Nx,c 
Roe ~ 10 5 yr! MWGal ! | ———“—\ || —**_ 
E ? 150MWGal! / | 1055 pc? 
Nan. ( o | mpg + m, | (rp/rx (11) 
100 /A10kms-! 11 Mo 1.0 
This equation implies that PTDEs (with larger rp/ri) are 
more frequent than FTDEs. More specifically, from our simu- 
lations of TDEs of 1 Mo stars by 10 Mo, encounters with b ~ 1— 
2 would result in eccentric bound orbits, which can circularize 
over time through tides. A more detailed rate calculation would 


involve integrating Equation 11 over a range of star and black 
hole masses. 
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Fig. 9: Post-disruption orbital eccentricities, e, for a star of mass 1 Mo (left) and 0.5 Mo (right) initially in a parabolic orbit, due to 
a BH of mass 10 Me (top) and 40 Mo (bottom), as a function of impact parameter b and MS age (H, abundance). The dashed line 
represents parabolic orbits and separates the parameter space into bound (eccentric) and unbound (hyperbolic) orbits. The reasoning 


for the trends is the same as that for Figure 8. 


The above rates are derived for uTDEs produced in dense 
stellar systems such as globular clusters and nuclear star clusters. 
However, other channels also contribute to  TDEs (see discus- 
sion in Perets et al. 2016), such as ultra-wide binaries and triples 
in the field (Michaely & Perets 2016, 2020) and post-natal-kicks 
leading to close encounter between a newly formed BH and a 
stellar companion (Michaely et al. 2016). These channels may 
give rise uTDEs in young environments and in the field. 


4.6. Summary of fits 


Below we present the post-disruption fits detailed in previous 
sections for use in globular cluster codes. m,4-,pH refers to the 
gas mass bound to the black hole after disruption, including the 
accreted mass. 


4.6.1. Masses of remnant star and BH 


(==) = 10min {-2.10(b-1.95p ee) > 0} (12) 

(=) E (==) (13) 
" * 

4.6.2. Spin of remnant star 

| T =) = 1 (-1-50b+52.00) (z) (1 _ Ar | (14) 


4.6.3. Validity of fits 


It should be noted that the density concentration factor, conc, 
only captures the full disruption of a star. Similarly defined fac- 
tors for each radius r, i.e. (p(r) /[p)!! 3 are required to calculate the 
distance at which mass beyond r is lost (see Ryu et al. 2020b). 
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Stars of higher masses tend to have much steeper density pro- 
files when compared to 1 Mo stars, thus Equation 4, which fits 
the mass loss, is not universal. In particular, the mass loss fit for 
2 Mo stars would involve an exponent of ~ 3, unlike our derived 
exponent of ~ 2.1 in Equation 4. Similar arguments follow for 
the fits of spins and orbits. 


5. Implications 


In this section, we discuss the general observational implications 
and the caveats of our study. uTDEs can potentially produce both 
unique and peculiar transient events as well as longer-lived sys- 
tems. 


5.1. Low-mass X-ray binaries 


Low-mass X-ray binaries (LMXBs) are interacting binary sys- 
tems which consist of compact objects (neutrons stars or black 
holes) accreting mass from their low-mass (x 1 Mo) stellar com- 
panions, thereby producing X-rays. A few hundred of LMXBs 
have been detected in the Milky Way Galaxy (e.g., Liu et al. 
2007; Avakyan et al. 2023), highlighting the prevalence of such 
systems. However, theoretical (e.g., Portegies Zwart et al. 1997; 
Kalogera 1999) and population synthesis (e.g., Podsiadlowski 
et al. 2003; Kiel & Hurley 2006) studies have not been able to 
match the rates of observed LMBXs through isolated binary evo- 
lution alone. The primary reason is that a low-mass companion 
has insufficient energy to expel a common envelope during the 
giant phase of the primary star, thereby resulting in the spiraling 
in of the companion and an eventual merger. Thus, a dynamical 
channel formation channel for LMXBs, where a low-mass star 
is tidally captured by a black hole in a globular/nuclear cluster, 
holds promise (e.g., Voss & Gilfanov 2007; Michaely & Perets 
2016; Klencki et al. 2017). 

Some of our simulations of encounters of low-mass stars 
with black holes are not close enough for significant mass loss 
and can result in highly eccentric bound orbits. For example, 
in the parabolic encounter of a 1 Mo star and a 10 Mo BH, for 
b ~ 1.0-1.5 (see Section 4), the fractional mass loss is ~ 1— 
10%, and the post-disruption bound orbit has a period of 0.1—1 yr 
and an eccentricity of ~ 0.97-0.98. However, tidal dissipation at 
subsequent periapsis passages can reduce the eccentricity and 
orbital period (e.g., Klencki et al. 2017) such that the resulting 
bound system evolves into a compact LMXB. 

Another type of system that is potentially formed through 
dynamical encounters are non-interacting BH-star systems, e.g., 
El-Badry et al. (2023a,b); Chakrabarti et al. (2023). These sys- 
tems consist of ~ 1 Mo stars orbiting ~ 10 Mo black holes in 
moderately-eccentric orbits with periods of ~ 1 yr. These param- 
eters correspond quite well with our simulations with b ~ 1.0— 
1.5, which lends credence to dynamical formation. 


5.2. Intermediate-mass black hole growth 


Intermediate-mass black holes (IMBHs; see Greene et al. 2020 
for a review), with masses greater than those of stellar-mass 
black holes (= 10? Mo) and less than those of supermassive 
black holes (z 109 Mo), are expected to exist in the centers 
of dwarf galaxies (e.g., Kunth et al. 1987; Filippenko & Sar- 
gent 1989; Reines 2022; Gültekin et al. 2022) and globular clus- 
ters (e.g., Davis et al. 2011; Farrell et al. 2014; Pechetti et al. 
2022). Among several IMBH formation mechanisms (see Volon- 
teri et al. 2021 for review), one possible avenue is runaway mass 
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growth through tidal capture and disruption, which has been ex- 
amined analytically (e.g., Stone et al. 2017) and numerically 
(e.g., Rizzuto et al. 2023; Arca Sedda et al. 2023). One of the 
most critical factors in this mechanism for BH growth is the 
amount of debris mass that ultimately accretes onto the BH. For 
example, Rizzuto et al. (2023) assume that a star approaching a 
BH within the nominal tidal radius r, is fully destroyed and 5096 
of the stellar mass is accreted onto the BH. 


While these assumptions may be valid for a part of the pa- 
rameter space, it is important to develop a prescription for the 
outcomes of TDEs that works over a broader parameter range 
to more accurately assess the possibility of BH mass growth 
through TDEs. Our simulations can provide such prescriptions 
that can improve the treatment of TDEs. Firstly, our simulations 
show, like previous work on TDEs by supermassive black holes 
(e.g., Guillochon & Ramirez-Ruiz 2013; Mainetti et al. 2017; 
Goicovic et al. 2019; Ryu et al. 2020a,b,c,d; Law-Smith et al. 
2020), that the periapis distances at which stars are fully or par- 
tially destroyed depend on the stellar internal structure. Ryu et al. 
(2020b) analytically demonstrated that the nominal tidal radius 
r, is a more relevant quantity for partial disruption events in- 
volving stars with m, z 1 Mo. Only for lower-mass stars, where 
Pconc ~ 1, does r; come close to the genuine full disruption ra- 
dius^. Our fitting formula for the fractional mass loss, Equation 
4, can provide a better prescription for determining the fate of the 
star (full or partial disruption) as well as the mass of the remnant. 


The amount of debris accreted onto the BH remains highly 
uncertain. To the zeroth order, if the accretion rate is super- 
Eddington, the strong radiation pressure gradient would drive 
strong outflows, hindering continuous and steady accretion 
(Sadowski et al. 2014). However, the accretion efficiency would 
be collectively affected by many factors, including magnetic 
fields, black hole spins, accretion flow structure, and jet for- 
mation (e.g., Sadowski et al. 2014; Jiang et al. 2019; Curd & 
Narayan 2023; Kaaz et al. 2023). Our fitting formulae for the 
fractional mass loss cannot provide an accurate prescription for 
the accreted mass but can place constraints on the maximum 
mass that can be accreted onto the BH. 


5.8. Fast blue optical transients 


Fast blue optical transients (FBOTs) are a class of optical tran- 
sients characterized by high peak luminosities > 10? ergs", 
rapid rise and decay times of the order of a few days, blue 
colors (Perets et al. 2011; Drout et al. 2014), and peak black- 
body temperatures of a few 10* K. Although a majority of these 
events can be explained by supernovae with low-mass ejecta 
(Pursiainen et al. 2018), a ‘luminous’ subset, e.g., AT2018cow 
(‘the Cow’; Smartt et al. 2018; Prentice et al. 2018), AT2018lIqh 
(Ofek et al. 2021), AT2020mrf (Yao et al. 2022), CSS161010 
(Coppejans et al. 2020), ZTF18abvkwla ('the Koala'; Ho et al. 
2020), AT2020xnd (‘the Camel’; Perley et al. 2021), AT2022tsd 
(‘the Tasmanian Devil’; Matthews et al. 2023), AT2023fhn (‘the 
Finch'; Chrimes et al. 2024), are too bright at their peaks and 
fade too rapidly to be explained by supernovae. Proposed mech- 
anisms involve compact objects, such as black hole accretion or 
magnetars formed in core-collapse supernovae (Prentice et al. 
2018), mergers between Wolf-Rayet stars and compact objects 
accompanied by hyper-accretion (Metzger 2022), or tidal dis- 


^ The genuine full disruption radius and the partial disruption radius 
are proportional to pe! ? and p 2 respectively (see Section 4 of (Ryu 
et al. 2020b)). 
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ruption events by intermediate-mass BHs (Kuin et al. 2019) or 
stellar-mass BHs (Kremer et al. 2021). 

Some luminous FBOTs also emit in radio or X-ray (e.g., 
Margutti et al. 2019; Ho et al. 2020) which indicates the pres- 
ence of a circumstellar medium that may result from partial 
TDEs (Kremer et al. 2021). Our simulations support this notion 
by showing the presence of debris near the BH with occasional 
weak outflows from the BH. In addition, the remnants in our sim- 
ulations are often bound to the BH and can subsequently undergo 
multiple partial disruptions, e.g., our simulations of TAMS 1 Mo 
stars and 10 Mo BHs for b = 0.25 or b = 0.33. While we cannot 
confirm a direct correlation between FBOTs and uTDEs, multi- 
ple disruption events of a star may fulfill some of the astrophys- 
ical conditions (e.g., the presence of a gas medium and rapid ac- 
cretion onto a compact object) necessary to explain FBOTs with 
radio emission. For such cases, the detection of repeated bursts 
in FBOTs could provide valuable constraints on their formation 
mechanisms. 


5.4. Ultra-long gamma-ray bursts and X-ray flares 


The exact appearance of transients resulting from pTDEs is 
highly uncertain and depends on various assumptions. Perets 
et al. (2016) suggested that, if the accretion onto the BH is ef- 
ficient and gives rise to jets, there is a possibility of producing 
ultra-long gamma-ray bursts (GRBs) and/or X-ray flares (XRFs). 

Assuming a proportional relation between the material ac- 
creted to the compact object and the luminosity produced by 
the jet, we can generally divide wTDEs into two regimes: (1) 
tmin >> fac; and (2) tmin < tacc, Where tmin is the typical fallback 
time of the debris following the disruption, and face is the typical 
viscous time of the accretion disk that can form around the BH 
(see Perets et al. 2016 for more details). 


— [gin >> tacc: When the mass of the star is much smaller than 
that of the BH, the accretion evolution is dominated by the 
fallback rate, i.e., the light curve should generally follow the 
regular power-law (e.g. t^?/? power-law for a full disruption). 

— tmin X tacc: When the masses of the star and the BH are 
comparable, the fallback material is expected to accumulate 
and form a disk on the fallback time, which then drains on 
the longer viscous time maintaining a low accretion rate. We 
would expect the flaring to begin only once the material is 
accreted onto the compact object. Therefore, we expect four 
stages in the light curve evolution: (1) A fast rise of the ac- 
cretion flare once the disk material is processed and evolves 
to accrete on the compact object. (2) Accretion from the disk 
until the accumulated early fallback material is drained; if we 
assume a steady state accretion until drainage, one might ex- 
pect a relatively flat light curve. (3) Once the disk drains the 
accumulated early fall-back material the light curve should 
drop steeply. (4) The continuous fall-back of material would 
govern the accretion rate at times longer than the viscous 
time, and the accretion rate should then follow the 15⁄3 rate 
(or a different power-law). The exact light curve at the early 
stages is difficult to predict, but we do note a non-trivial sig- 
nature of the wTDEs relating to the early and late stages. 


The expected properties of wTDEs are consistent with and 
might explain the origins of ultra-long GRBs (Levan et al. 2014): 
long-lived (^ 10^ s) GRBs which show an initial plateau fol- 
lowed by a rapid decay. uTDEs may also explain the origin of 
some Swift TDE candidates (Bloom et al. 2011; Burrows et al. 
2011; Cenko et al. 2012; Brown et al. 2015), suggested to be 


produced through a TDE by a supermassive BH. The typical 
timescales for the latter are longer than the observed 10° s, chal- 
lenging the currently suggested origin, but quite consistent with 
a UTDE scenario. 

HTDEs producing ultra-long GRBs are also expected to pro- 
duce afterglow emission if/when the strong outflows and jets 
launched by accretion propagate into the surrounding medium 
and the resulting shock interaction produces high energy parti- 
cles which give rise to synchrotron and inverse Compton radi- 
ation. Such observational signature of uTDE has been little ex- 
plored, though analogous observation and modeling have been 
suggested for TDEs by massive black holes (e.g., Giannios & 
Metzger 2011; Zauderer et al. 2011). Detection of afterglows 
from uTDEs would provide information on the energetics and 
dynamics of the outflow, robust independent evidence for the jet 
collimation, as well as identification of the host galaxy. 


5.5. Remnants 


Partial uL TDEs of stars by stellar BHs can spin up the stars sig- 
nificantly due to strong tidal interactions (see Fig.7 and Alexan- 
der & Kumar 2001). Though magnetic breaking can potentially 
spin down such stars over time, observations of highly spun-up 
low-mass old stars in globular clusters could provide potential 
signatures for partial uTDEs. In addition, the remnant undergoes 
violent chemical mixing during the first pericenter passage and 
has higher entropy than an ordinary star of the same mass and 
age (Ryu et al. 2020c, 2023a). If the mass loss is significant, a 
unique chemical composition profile may persist even after the 
remnant has relaxed to a stable state. This suggests that the rem- 
nant could exhibit unique asteroseismic signatures that distin- 
guish it from ordinary stars (Bellinger et al. 2023, in prep). 


5.6. Caveats 


As mentioned in Section 2, we exclude the effects of relativity 
and magnetic fields in our study. The Newtonian approximation 
is justified for TDEs by SBHs because most of the debris stays 
outside rp, which is approximately four orders of magnitude 
greater than the gravitational radius rs = Gmpy /c? >. However, 
relativistic effects and proper treatment of radiation are more im- 
portant for the accretion flow, which is beyond the scope of our 
paper. The effect of magnetic fields on the dynamics of the rem- 
nant is unclear, which we will leave for our future work. A final 
point to note is that, unlike globular clusters, nuclear star clus- 
ters have typical velocity dispersions of © = 100kms Figer 
et al. 2003; Schódel et al. 2009) depending on the galaxy, and 
subsequently, |1 — e| = 1072. Thus, encounters are not always 
parabolic, and for a more thorough study, hyperbolic encounters 
must be taken into account. 


6. Summary and conclusion 


We performed a grid of 58 hydrodynamics simulations, using the 
moving-mesh code AREPO, of partial (PTDE) and full (FTDE) 
tidal disruptions of main-sequence (MS) stars with stellar-mass 
black holes (SBH) on initially-parabolic orbits. Our varied pa- 
rameters include stellar mass m, (0.5 Mo and 1.0 Mo), SBH 
mass (10 Mo and 40 Mo), core hydrogen abundance He (a proxy 
for MS age f.) and impact parameter b = rj/ri. Our stellar mod- 
els are initialized in accordance with 1D detailed non-rotating 


5 On the other hand, for TDEs with SMBHs, r, and r, can be compa- 
rable in scale. 
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stellar profiles from MESA. We also define a star’s density con- 
centration parameter Pconc = (pc Ip)" 5. Our main results are sum- 
marized below: 


— The mass of a post-disruption stellar remnant m, decreases 
with decreasing b. A higher Peonc (more centrally concen- 
trated) results in a lesser mass loss for the same b. m, de- 
pends only on b and pconc, and not on mpy. Roughly half of 
the disrupted mass remains bound to the black hole, while 
the other half is unbound from the system. 

— The spin angular momentum of a stellar remnant L, in- 
creases after periapsis passage. For large b (when the mass 
loss is small), L, increases with decreasing b, but this trend 
reverses for very small b (close to FTDE) due to significant 
mass loss. This reversal occurs for lower b when conc is 
higher. The final spin is also independent of mgy. For low 
b values, the rotational velocities can be very close to break- 
up velocities. 

— Unlike mass and spin, the orbit of the remnant star also de- 
pends on mpu. For large b, eccentricity e and semimajor axis 
a decrease with decreasing b keeping a constant periapsis 
distance, i.e., the remnant tends to be more ‘bound’. Depend- 
ing on mpy and peonc, this no longer holds for approaches 
close to FTDE - e can increase and the remnant can become 
unbound. In general, a higher mgy/m, ratio and a lower oconc 
tend to be unbinding. 

— We provide relatively simple and accurate fitting formulae 
for the masses and spins of the remnant stars. These analyti- 
cal fits are listed in Section 4.6 for easy access. 


We discussed the implications that u TDEs can have on low- 
mass X-ray binary formation, intermediate-mass black hole for- 
mation, and transient searches. Given the limited range of stellar 
masses and ages that we covered in this work, we will extend 
our suite of simulations in the future to cover a wider parame- 
ter space, including massive stars. It is necessary to explore this 
regime to model stellar dynamics in dense stellar clusters. 
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Appendix A: Summary of simulation results 


Table A.1: Post-disruption parameter values for our 58 simulations. 


My MBH b He Ly / 10% Mo Tl, BH Munbound Lo / 1045 a e 
[Mo] [Mo] [gcm°s!] [Me] [Me] [Ms] [gem?s] [au] 
1 0.5 10.0 0.70 0.25 2 0.00 0.29 0.21 N/A N/A N/A 
2 0.5 10.0 0.70 0.33 2 0.00 0.22 0.28 N/A N/A N/A 
3 0.5 10.0 0.70 0.50 2 0.00 0.23 0.17 N/A N/A N/A 
4 0.5 10.0 0.70 0.75 2 0.15 0.21 0.14 87 -0.07 1.06 
5 0.5 10.0 0.70 1.00 2 0.25 0.14 0.11 97 -0.27 1.02 
6 0.5 10.0 0.70 1.50 2 0.46 0.03 0.01 117 1.91 0.995 
7 0.5 10.0 0.70 2.00 2 0.50 0.00 0.00 17 2.61 0.996 
8 0.5 10.0 0.70 2.50 2 0.50 0.00 0.00 1 42.3 1.00 
9 0.5 40.0 0.70 0.25 2 0.00 0.27 0.23 N/A N/A N/A 
10 0.5 40.0 0.70 0.33 2 0.00 0.30 0.20 N/A N/A N/A 
11 0.5 40.0 0.70 0.50 2 0.00 0.30 0.20 N/A N/A N/A 
12 0.5 40.0 070 0.75 2 0.00 0.31 0.19 N/A N/A N/A 
13 0.5 40.0 0.70 1.00 2 0.21 0.18 0.11 76 -0.35 1.03 
14 0.5 40.0 0.70 1.50 2 0.46 0.03 0.01 110 -0.11 1.00 
15 0.5 40.0 0.70 2.00 2 0.50 0.00 0.00 15 11.4 0.998 
16 0.5 40.0 0.70 2.50 2 0.50 0.00 0.00 1 161 1.00 
17 1.0 10.0 0.00 0.25 12 0.04 0.42 0.54 2 0.02 0.916 
18 1.0 10.0 0.00 0.33 12 0.33 0.27 0.40 101 0.04 0.923 
19 1.0 10.0 0.00 0.50 12 0.68 0.11 0.21 410 0.10 0.945 
20 1.0 10.0 0.00 0.75 12 0.89 0.05 0.06 354 0.32 0.972 
21 1.0 10.0 0.00 1.00 12 0.97 0.02 0.01 199 1.15 0.989 
22 1.0 10.0 0.00 1.50 12 1.00 0.00 0.00 49 13.2 0.999 
23 1.0 10.0 0.00 2.00 12 1.00 0.00 0.00 10 54.6 1.00 
24 1.0 10.0 0.34 0.25 10 0.00 0.47 0.53 N/A N/A N/A 
25 1.0 10.0 0.34 0.33 10 0.09 0.47 0.44 99 -0.06 1.12 
26 1.0 10.0 0.34 0.50 10 0.49 0.20 0.31 245 0.08 0.950 
27 1.0 10.0 0.34 0.75 10 0.83 0.09 0.08 440 0.21 0.965 
28 1.0 10.0 0.34 1.00 10 0.95 0.03 0.02 250 0.64 0.984 
29 1.0 10.0 0.34 1.50 10 1.00 0.00 0.00 49 7.77 0.998 
30 1.0 10.0 0.34 2.00 10 1.00 0.00 0.00 10 87.3 1.00 
31 1.0 10.0 0.70 0.25 16 0.00 0.37 0.63 N/A N/A N/A 
32 1.0 10.0 0.70 0.33 16 0.00 0.49 0.51 N/A N/A N/A 
33 1.0 10.0 0.70 0.50 16 0.34 0.29 0.37 163 0.28 0.986 
34 1.0 10.0 0.70 0.75 16 0.68 0.16 0.16 440 0.2.2 0.968 
35 1.0 10.0 0.70 1.00 16 0.90 0.06 0.04 387 0.35 0.974 
36 1.0 10.0 0.70 1.50 16 1.00 0.00 0.00 67 3.43 0.996 
37 1.0 10.0 0.70 2.00 16 1.00 0.00 0.00 8 25.5 0.999 
38 1.0 40.0 0.00 0.25 12 0.07 0.45 0.48 56 -0.05 1.14 
39 1.0 40.0 0.00 0.33 12 0.28 0.27 0.45 66 0.17 0.969 
40 1.0 40.0 0.00 0.50 12 0.63 0.16 0.21 340 0.58 0.984 
41 1.0 40.0 0.00 0.75 12 0.88 0.06 0.06 332 2.24 0.994 
42 1.0 40.0 0.00 1.00 12 0.96 0.02 0.02 192 9.63 0.998 
43 1.0 40.0 0.00 1.50 12 1.00 0.00 0.00 50 102 1.00 
44 1.0 40.0 0.00 2.00 12 1.00 0.00 0.00 8 153 1.00 
45 1.0 40.0 0.34 0.25 10 0.00 0.42 0.58 N/A N/A N/A 
46 1.0 40.0 0.34 0.33 10 0.00 0.39 0.61 N/A N/A N/A 
47 1.0 40.0 0.34 0.50 10 0.42 0.22 0.36 169 0.95 0.992 
48 1.0 40.0 0.34 0.75 10 0.81 0.09 0.10 413 1.84 0.993 
49 1.0 40.0 0.34 1.00 10 0.95 0.03 0.02 245 5.09 0.997 
50 1.0 40.0 0.34 1.50 10 1.00 0.00 0.00 50 62.6 1.00 
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51 1.0 40.0 0.34 2.00 10 1.00 0.00 0.00 10 180 1.00 
52 1.0 40.0 0.70 025 16 0.00 0.55 0.45 N/A N/A N/A 
53 1.0 40.0 0.70 033 16 0.00 0.50 0.50 N/A N/A N/A 
54 1.0 40.0 0.70 0.50 16 0.21 0.39 0.40 156 -0.12 1.06 
55 1.0 40.0 0.70 0.75 16 0.68 0.18 0.19 365 -4.17 1.00 
56 1.0 40.0 0.70 1.00 16 0.89 0.06 0.05 368 2.9] 0.995 
57 1.0 40.0 0.70 1.50 16 1.00 0.00 0.00 66 29.4 0.999 
58 1.0 40.0 0.70 2.00 16 1.00 0.00 0.00 8 119 1.00 


Notes. The columns, respectively, denote star mass m,, BH mass mpg, impact parameter b, central hydrogen abundance H., initial (almost non- 
rotating) star spin angular momentum L,, remnant mass m,, gas mass bound to BH m,_,gy, gas mass unbound from the BH-star system bouna. 


remnant spin angular momentum L,, orbital semimajor axis a and orbital eccentricity e. 
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